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C\J ABSTRACT 

> ; 

0^ , We present submillimetre photometry of 11 SCR and 6CE radio quasars of similar 

■ 151-MHz radio luminosity and redshifts to the radio galaxies studied at 1.3 < z < 2 
I by Archibald et al. (2001). We detect all bar one quasar at 850 (a.m, and five quasars 

^ . are confirmed as dusty hyperluminous objects by detections at shorter wavelengths. 

I ' Our observations reveal a clear difference between the submillimetre luminosity dis- 

\ tributions of the radio quasars and a matched sample of radio galaxies: the quasars 

■ are ~ 5-times brighter than the radio galaxies and ^ 2-times brighter accounting for 
' possible synchrotron contamination. This difference is in quantitative agreement with 

a receding-torus unified scheme for radio sources in which the torus opening angle de- 
Q ■ pends on quasar optical luminosity, provided that there is a close relationship between 

I the optical and submillimetre luminosities. The implication is that quasar-heated dust 

C/3 , dominates the submillimetre emission for powerful radio quasars at z ~ 1.5. We find 

' a significant anti-correlation between submillimetre/far-infrared luminosity LpiR and 

^ \ radio source age in the sense that hyperluminous quasars tend to be associated with 

•'H , young (< 10^ yr) sources. This suggests that the processes controlling LpiR are syn- 

' chronized with the jet-triggering event and that ipiR declines on a ~ 10^ yr timescale. 

^ , There is evidence for a small amount of obscuration in the hyperluminous quasars from 

■ - - i reddening and associated or broad absorption lines. We speculate that shortly after 

jet triggering, dust is more widely distributed than at later times and is heated by the 
quasar nucleus and/or a short-lived synchronized starburst. Any such starburst would 
form only a few per cent of the total stellar mass in agreement with the evidence for 
dominant old stellar populations in z ^ 1.5 radio galaxies. In contrast, high-redshift 
(z > 3) radio galaxies with similar submillimetre luminosities could have longer star- 
formation timescales due to the greater availability of gas and could be forming the 
bulk of their stellar populations. 

Key words: galaxies: active - galaxies: evolution - quasars: individual: 3C 268.4, 
3C 298, 3C 318, 6C 1045+3513 



1 INTRODUCTION 

Submillimetre observations provide an excellent method of 
detecting dusty star-forming galaxies at high redshifts: UV 
radiation from hot stars is absorbed by dust, re-radiated 
in the far-infrared and redshifted into the submillimetre. 
The steep submillimetre spectral index leads to large k- 
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corrections which, at high redshifts, compensates for increas- 
ing cosmological dimming so that an object with a given sub- 
millimetre luminosity has a similar observed flux-density at 
z = 1 to a similarly-luminous one at 2 = 5 (Blain & Longair 
1993). 

Powerful radio sources can be easily identified out to 
high redshifts and are thought to be associated with mas- 
sive elliptical galaxies, or their progenitors. The existence of 
this close association is well known at low-to-intermediate 
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redshifts (e.g. Owen & Laing 1989) and it is likely to be 
maintained to higher redshifts by scaling relations between 
radio luminosity, accretion luminosity, black hole mass and 
(dark matter) halo mass (e.g. Rawlings & Saunders 1991; 
Magorrian et al. 1998). At low-redshifts, radio galaxies ap- 
pear to be kinematically-relaxed massive ellipticals (McLure 
et al. 1999; Dunlop et al. 2001; Bettoni et al. 2001) with old 
stellar populations (Nolan et al. 2001). These facts, together 
with the small scatter in the near-infrared Hubble diagram 
to z ~ 3 (Jarvis et al. 2001a), suggest that most of the 
star formation in powerful radio galaxies occurred at early 
cosmic epochs. 

In the first large study of the submillimetre emission of 
high-redshift radio galaxies, Archibald et al. (2001; hereafter 
AOl) made sensitive (850 \i.m sensitivity ~ 1 mjy rms) sub- 
millimetre observations of 47 radio galaxies with the SCUBA 
array on the JCMT. Their sample covered the redshift range 
1 < z < 4.5 and they detected dust emission in 14 objects. 
The main result of their study was a clear difference between 
the detection rate at 2: < 2.5 of 15% and the detection rate 
at z > 2.5 of 75%. This indicates that the submillimetre 
luminosity of radio galaxies increases systematically with 
redshift, suggesting that the bulk of star formation in radio 
galaxies, and by inference in all massive galaxies, occurred 
at epochs corresponding to high redshifts. 

Unified schemes for radio galaxies and radio 
quasars|^ropose they are the same objects viewed along 
different lines-of-sight, such that in radio galaxies the nu- 
clear regions are obscured by a dusty torus or warped disc 
(Scheuer 1987; Barthel 1989). Such theories have proved 
remarkably successful (see Antonucci 1993 for a review), at 
least at high radio luminosities (Willott et al. 2000a). For 
objects at redshifts 2 ~ 1.5, the observed 850 \i.m radiation 
was emitted at a rest-frame wavelength of « 300 \im. At 
these wavelengths, dust is expected to be optically thin 
and its thermal emission isotropic (e.g. Granato & Danese 
1994; Efstathiou & Rowan- Robinson 1995). Therefore, ac- 
cording to the simplest type of radio-loud AGN unification 
scheme, the submillimetre properties of radio galaxies and 
radio quasars, matched in terms of low- frequency radio 
luminosity and redshift, should be similar. In more realistic 
unified schemes, however, they may not be expected to 
be identical. Rawlings & Saunders (1991) and Simpson 
(1998) have argued that samples of radio galaxies and radio 
quasars, matched in one orientation-independent quantity 
(i.e. low- frequency radio luminosity), are expected to have 
different distributions in a second orientation-independent 
quantity (i.e. [O III] emission line luminosity) if the two 
quantities are linked by a scaling relation with intrinsic 
scatter. The first motivation behind the work described 
in this paper was to explore unified schemes further by 
looking for significant differences between the submillimetre 
properties of radio galaxies and radio quasars matched in 
151-MHz luminosity and redshift. 

The second motivation behind this paper was to ex- 

^ We will use the term 'radio quasars' to refer to lobe-dominated 
steep-spectrum radio sources associated with quasar nuclei. These 
are the dominant population of radio-loud quasars in low- 
frequency-selected samples such as 3CRR and 6CE, whereas com- 
pact flat-spectrum quasars dominate in high-frequency-selected 
samples. Flat-spectrum quasars often have high submillimetre 



plore further the suggestion of Blundell & Rawlings (1999) 
that some part of the evolution in the submillimetre proper- 
ties of radio galaxies with redshift is a selection effect linked 
to the inevitable youth of the highest-redshift radio sources 
from flux-limited samples. AOl have already considered this 
suggestion, arguing that it is not a primary driver of the 
observed redshift dependencies. However, with a small dy- 
namic range in source age, and somewhat uncertain selec- 
tion effects, the high-redshift sample of Archibald et al. is 
not ideally suited to a test of the idea that submillimetre 
properties might be correlated with the time since the jet- 
triggering event. 

At z ^ 1.5 AOl have already obtained submillimetre 
data on an unbiased sample of 3CRR and 6CE galaxies, 
and we aimed to complement these with new observations 
of a matched sample of radio quasars. In Section 2 we de- 
scribe the observations and results. In Section 3 we consider 
the likely contribution to the measured submillimetre fluxes 
from an extension of the radio synchrotron spectrum. In Sec- 
tion 4 we consider explanations for our results. We present 
our conclusions in Section 5. The convention for all spectral 
indices a is that flux-density Sv oc v~°' . We assume through- 
out that Ho = 70 km s^^Mpc"\ Q,m = 0.3 and = 0.7. 



2 OBSERVATIONS 

We have defined a sample of radio quasars in the redshift 
range 1.37 < z < 2.0 selected from two complete radio sam- 
ples (3CRR - Laing, Riley & Longair 1983; 6CE - Rawlings, 
Eales & Lacy 2001) to compare with the radio galaxies from 
the same samples in the same redshift range in the study of 
AOl. Due to scheduling reasons only 3CRR quasars within 
the Right Ascension range 04'' < RA < 22'' are included 
in our sample. These were supplemented with two quasars 
from the revised 3C catalogue (3CR - Spinrad et al. 1985). 
These two quasars, 3C 298 and 3G 280.1, are not included 
in the 3CRR sample due to reasons of low declination and 
radio flux-density just below the 3CRR limit, respectively. 
They are the only high galactic latitude quasars in the target 
redshift range in the 3CR sample which are not included in 
the 3CRR sample. The total sample consists of 11 quasars 
and the comparison sample extracted from AOl contains 12 
radio galaxies (Table 1). 

The sample of radio quasars was observed in photom- 
etry mode with the SCUBA bolometer array at the JCMT 
in March, April, October 1999, November 2000, February 
and April 2001. Observations were made simultaneously at 
850 |i.m and 450 \uxi and the sensitivities reached at 850 
\im are ^ 1 mJy rms. The atmospheric zenith opacity at 
850 \x.m ranged from 0.16 to 0.52 with a median of 0.22. 
The data were reduced using the SURF package, following 
the same method as in AOl for consistency. The resulting 
flux-densities are given in Table 1 along with the data for 
the matched radio galaxy sample. Ten out of eleven radio 
quasars were detected at 850 \un at > 2a significance (eight 
at > 3cr), compared with only four > 2(t detections out of 
the twelve radio galaxies. At the shorter wavelength of 450 
\im, three quasars have > 2a detections, compared to none 

fluxes simply by virtue of a continuation of their synchrotron 
spectra from radio to shorter wavelengths. 
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13.6 + 18.0 
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-19.3+ 8.2 
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6C 0905+3955 
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934.5 


3.62 + 0.89 
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27.53 


6.0 
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1.779 


27.61 


1.0 


36.45 


435.2 


0.16 + 1.25 


45.1 + 26.6 


2.66 


98 



Table 1. Submillimetre observations and basic data for the sample of radio quasars (class 'Q') and radio galaxies ('G') in the redshift 
range 1.37 < z < 2.0 described in Sec. 2. Detections at the ^ 2a level are given in bold type. Submillimetre flux-deusities of the radio 
galaxies are from AOl. No corrections have been made for possible contributions of synchrotron emission to the submillimetre fluxes. 
Radio core fluxes are quoted at 5 GHz (observed-frame) in most cases. Radio data for the radio galaxies and quasars are from the 
references listed in AOl and Blundell et al. (in prep.) respectively. Data for the [OIII] emission line luminosities have been taken from 
Jackson & Rawlings (1997) for 3CR sources and Rawlings et al. (2001) for 6C sources, apart from these exceptions: 3C 241 - Hirst, 
Jackson & Rawlings (2002); 3G 270.1 and 6G 0955+3844 - Willott et al. (in prep.); 3G 298 - Espey et al. (1989); 3C 318 - Willott et al. 
(2000b); 6C 1045+3513 (this paper). For some sources, [OIII] has not been observed and a different narrow line flux is used and the [OIII] 
flux estimated from the Une ratios in McCarthy (1993). These sources are 3C 298 ([OH]): 6C 0820+3642 (Lya); 6C 0901+3551(Lya); 
60 0905+3955 (Lya); 6C 0919+3806 (CII]); 6C 1045+3513 ([NeV]); 60 1204+3708 (Lya). For the quasar 30 432, no narrow line data 
have been published and the [OIII] flux has been estimated assuming a rest- frame equivalent width of 30 A (e.g. Miller et al. 1992). 



of the radio galaxies (although the quasar observations are 
generally more sensitive at 450 than those of the radio 

galaxies due to the introduction of new filters in October 
1999). The mean quasar SsrM is 6.3 ± 1.0 mjy and the moan 
radio galaxy Ss50 is 1.2 + 0.6 mJy. Therefore the quasars are 
brighter than the radio galaocies in the submillimetre by a 
factor of 5. 

Relating the observed submillimetre flux-density to a 
submillimetre/far-infrared luminosity is not a trivial matter, 
given the poor constraints on the shape of the dust spectra 
for most sources which have only been detected at one or two 
wavelengths. Therefore, we adopt a template dust spectrum 
as has been determined from observations of a sample of 
similarly luminous (submillimetre/far-infrared luminosities 
ipiR ~ lO^i©) objects by Priddey & McMahon (2001). 
This template consists of isothermal optically-thin dust with 
temperature T = 41 K and emissivity index l3 — 1.95. For 
our adopted cosmological model and dust template, the ob- 
served 850iJ.m flux-density Ss50 of a source with a given Lfir 
does not change appreciably over the redshift range z = 1 
to z = 5. The conversion between Lfir and Sgso for our 
adopted dust spectrum is Lfir ~ 2 x 10^^(5'85o/mJy)LQ. 
However, for higher dust temperatures or multi-temperature 
models (which are required to fit the infra-red SEDs of most 
luminous AGN - e.g. Rowan- Robinson 2000), Lfir could 



be up to an order of magnitude greater than this, so these 
values of Lfir should be treated as lower limits. Therefore 
the radio quasars we have detected with SCUBA have lower 
limits on their far-infrared luminosities in the range 4 x 10^^ 
to 2 X 10^^ Lq. It is likely that many of them would be clas- 
sified as hyperluminous infrared sources (Lfir ^ IO'^Lq) if 
their far- and mid-infrared fluxes were measured. Through- 
out most of this paper we will use the observed submillimetre 
flux-density in preference to the far-infrared luminosity and 
thereby avoid uncertainty in Lfir due to the unknown dust 
spectral shape. 

Five of the radio quasars are also detected at 
wavelengths shorter than 850n.m. 3C 191, 3C 318 and 
6C 1045+3513 are detected at 450^m with SCUBA and 
3C 268.4, 3C 298 and 3C 318 are detected in the far-infrared 
by ISO and/or IRAS. These objects are discussed in further 
detail here. 

2.1 Notes on objects detected at 100 — 450iJ.m 
3C 191 

The submillimetre spectral index of 3C 191 determined 
from the 850n.m and 450n.m fluxes is —2.4 ± 0.6, consistent 
within 2a with the value of —3.2, expected for our assumed 
dust template. We find that Lfir > 1.2 x lO" Lq. The op- 



4 



Willott et al. 



Rest Wavelength / /^m 



Rest Wavelength / /j^m 



10" 
10" 



3 10" 



X 10" 



10 

10" 




Rest Frequency / Hz 



16 



Figure 1. Rest-frame SED from radio to optical of the liyperlu- 
minous IR quasar 3C 268.4. For observations which separate the 
core emission from that of the extended radio components, the 
core flux is shown as an open circle. The dotted line is a a = 1 
core spectrum which dominates the millimetre flux. The dashed 
curve is the best fit thermal dust spectrum (T = 48 K, /3 = 2.5, 
optically-thin at A > 50i^m). The solid curve shows this added to 
the synchrotron core spectrum. 
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Figure 2. Rest-frame SED of the hyperluminous IR quasar 
3C 298 — symbols as for Fig. |^. The dotted line is a a = 1 core 
spectrum which dominates the millimetre flux. The solid curve 
shows this added to the thermal IR spectrum which is modelled 
as a combination of greybodies of temperatures ranging from 50 
K (this component is shown by a dashed line) to 1300 K. The 
dot-dashed line shows that the submillimetre/far-infrared data 
points could equally well be fitted by a T=100 K dust spectrum. 



tical spectral slope is Oopt ~ 0.7 (Barthel, Tytler & Thom- 
son 1990), quite similar to the typical radio quasar value of 
Oopt = 0.5 (Brotherton et al. 2001). This indicates that there 
is little reddening of the quasar light along our line-of-sight. 
However, 3C 191 has very strong associated absorption lines 
with a CIV absorption equivalent width of 6.1 A (Anderson 
et al. 1987). Hamann et al. (2001) studied these absorption 
systems with high-resolution spectroscopy and determined 
the distance of the absorbers from the nucleus to be ~ 30 
kpc with a corresponding characteristic flow time of 3 x 10^ 
years. The distances of the radio hotspots from the nucleus 
are of a similar size-scale (projected linear size=41 kpc), 
corresponding to an approximate age of the radio source of 
3 X 10^ years (see Sec. 4.2). 

3C 268.4 

3C 268.4 is detected at 850n.m (5.1±1.2 mjy) but not at 
450[i-m. Recent detections at 170^i.m and 90vim with ISO and 
at 1.3 mm with IRAM (6.1 ± 1.6 mJy) indicate a high far- 
IR luminosity and some synchrotron contamination of the 
mm fluxes (Andreani et al. 2001). Estimating the millimetre 
contamination by synchrotron emission is not a trivial mat- 
ter due to the lack of other data at frequencies higher than 
15 GHz (Fig. |l|). For reasons to be discussed in Sec. |^ we 
expect the extended radio lobe spectrum to steepen sharply 
at high frequencies and the main synchrotron component in 
the millimetre regime to be the compact core. It is also clear 
in Fig. |l| that the core spectrum also steepens between 15 
GHz and 1.3 mm. To estimate the synchrotron contamina- 
tion to the 850^tm flux-density we assume a millimetre core 
spectrum with a — 1 which gives Sgso = 3 mJy. As we will 
show below, the properties of the dust spectrum require a 



contribution to the 850(i.m flux-density from dust which is 
in approximate agreement with this contribution from the 
core. 

We have fitted the far-IR and sub-mm data with a 
range of thermal dust spectra. Using only the ISO fluxes 
and the 450ii.m upper limit, the strongly peaked spectrum 
at ~ 170fi.m (observed frame) provides strong constraints 
on the dust temperature and emissivity index (Fig. ^). The 
best fitting values are T = 48 K and /3 — 2.5 for dust which 
is optically-thin at A > 50(i.m and this spectrum gives an 
expected 850 jxm flux-density of 2.7 mJy. The far-IR lumi- 
nosity of this dust spectrum is 2 x 10^'' Lq. This value for the 
emissivity index is larger than usually observed (Priddey & 
McMahon 2001; Dunne et al. 2000) but requiring /3 to be 
no larger than 2.0 gives a similar far-IR luminosity and a 
higher expected 850ti.m fiux-density of 4.0 mJy. Therefore 
we conclude that the dust contribution to the 850p.m fiux- 
density is significant and independent of its actual value, 
Lfir > 2 X IO^'^Lq and the object is classified as hyperlu- 
minous. 3C 268.4 is not heavily reddened in the rest-frame 
UV but does show strong C IV associated absorption (equiv- 
alent width > 1.9 A; Anderson et al. 1987). 

3C 298 

The source with the highest 850 [im fiux-density in our 
sample is 3C 298. This quasar is extremely bright at all wave- 
lengths from the radio through to X-rays. It was detected 
by IRAS at 60 \Lm and by ISO at five wavelengths ranging 
from 12.8 to 170 ^m (Meisenheimer et al. 2001). 3C 298 has 
a strong radio core (515 mJy at 5 GHz; Akujor et al. 1991) 
but is lobe-dominated at GHz frequencies. The measured 
1.3 mm flux-density of 14.1 mJy (Meisenheimer et al. 2001) 
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shows that the core is still strong at millimetre wavelengths 
and makes a substantial contribution to the observed 850 
^.m flux-density. We have fitted the SED of 3C 298 to at- 
tempt to determine the relative contributions of dust and 
synchrotron emission (Figure |^ . This procedure is not triv- 
ial because the core spectrum appears to be fairly flat up 
to a rest-frame frequency of 200 GHz and then steepens 
sharply to a « 2 at higher frequencies. However, the flux- 
density measurements at these high frequencies have large 
errors and are only la from being consistent with a — 1, 
as is more common above the break in the core spectra of 
lobe-dominated quasars (e.g. Antonucci, Barvainis & AUoin 
1990; van Bemmel & Bertoldi 2001). Therefore we assume 
a core with a = 1 over the observed 1300 to 850 ^im range. 
We conclude that « 8 mjy of the 850 [im flux-density is due 
to synchrotron radiation and ~ 13 mJy is due to thermal 
dust emission. The reason why this quasar was not detected 
at 450 ^m was due to the poor sensitivity caused by the 
high airmass of the observation. The submillimetre and far- 
infrared data can be equally well fit by thermal spectra with 
temperatures T = 50 K (typical of starbursts) and T = 100 
K. The total IR luminosity of 3C 298 is 8 x lO^L©. 

3C 298 is a small radio source with projected linear size 
of 12.7 kpc. It shows very strong CIV associated absorption 
with equivalent width 4.5 A (Anderson et al. 1987). The 
optical spectral slope is Oopt = 1.5, which is redder than 
typical quasars. This reddening is most likely due to dust 
associated with the absorbing gas. 

3C 318 

In Willott, Rawlings & Jarvis (2000b) we have already 
studied the spectral energy distribution (SED) of 3C 318 and 
determined a total IR luminosity of > 7x IO^'^Lq (corrected 
for the cosmology used in this paper). The quasar optical 
spectrum is found to be reddened, indicating the presence 
of dust in the quasar environment. 3C 318 is a small (6.8 
kpc) radio source. 

6C 1045+3513 

6C 1045-1-3513 is detected at both 850 and 450 ^m, with 
a submillimetre spectral index of — 1.3±0.7. On the basis of 
the 850 \im observations we find Z/pir, > 2 X 1O^^L0. In the 
Appendix we present new optical spectroscopy of this ob- 
ject, showing it to be a reddened quasar, similar to 3C 318. 
The spectrum also shows a CIV broad absorption feature. 
6C 1045-1-3513 is the smallest radio source of the entire sam- 
ple of quasars and radio galaxies with a projected linear size 
of just 1.7 kpc. 



2.2 Comparison of radio quasars and radio 
galaxies 

Fig. 1^ plots the observed 850 \im flux-density S'gso against 
low-frequency (151 MHz) radio luminosity Lisi for the 
matched quasar and radio galaxy samples. The main point 
to note from Fig. |^ is that there is a clear difference in 
the flux-density distributions of the quasar and radio galaxy 
samples, with the mean quasar submillimetre flux brighter 
by a factor of 5j^ 

^ A second point to note from Fig. ^ is that there is only marginal 
evidence for a correlation between L151 and Sgso- The survival 
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Figure 3. 850 |^m flux-density Ssso against low-frequency (151 
MHz) radio luminosity L151 for the 1.37 < z < 2.0 radio quasars 
(filled circles) and radio galaxies (open circles). 2a upper limits 
are denoted by symbols with arrows. The five quasars with detec- 
tions at shorter wavelengths (100 — 450 [jm) have crosses through 
them. 3C 298 and 3C 268.4 are shown as two symbols joined 
by a line, indicating the total flux-density and that corrected for 
non-thermal emission. Also plotted are radio galaxies at higher 
redshifts (2 < 2 < 4.5) from AOl (small squares). 



Due to the large number of 850 \im upper limits for 
the radio galaxies, we use survival analysis statistical tests 
which can account for these limits (Feigelson & Nelson 1985; 
Isobe, Feigelson & Nelson 1986). Using a variety of tests to 
compare the 850 vim flux distributions of the quasars and 
radio galaxies (the Gehan, logrank and Peto-Prentice tests), 
the flux-density distributions are found to be different at 
> 99.9% signiflcance. Given the high fraction of upper limits 
for the radio galaxies, we cannot unequivocally determine 
the factor by which the quasars are brighter than the radio 
galaxies in the submillimetre. Using the above tests we flnd 
that the quasars are brighter than the radio galaxies by a 
factor 2 at the 95% confidence level. 



3 THE SYNCHROTRON CONTRIBUTION TO 
SUBMILLIMETRE FLUXES 

We have found the robust result that radio quasars have 
2 — 5-times higher submillimetre fluxes than radio galaxies 
of the same radio luminosity and redshift. If this emission is 
due to optically-thin dust, then this is inconsistent with the 
notion of radio quasars and radio galaxies being identical 
apart from some orientation dependence. Given that the ra- 
dio sources we are discussing are some of the most powerful 



analysis probability for a correlation being present is 80%. Re- 
moving the outlying point of 6C 1045+3513 increases the corre- 
lation probability to 95%. This matter has been discussed by 
Archibald et al. (2001), and we defer any further discussion to a 
wider investigation of the dependences of submillimetre luminos- 
ity on radio luminosity and redshift which will use new SCUBA 
data on radio sources from fainter samples than the SCR and 6CE 
samples used here (Rawlings et al., in prep.). 
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known, it is essential that we consider whether the observed 
submillimetre fluxes could be affected by an extrapolation 
of the nonthermal radio synchrotron emission. Radio lobes 
and hotspots generally have a very steep spectrum at mil- 
limetre wavelengths and are unlikely to contribute signif- 
icantly for many sources. This issue is discussed in detail 
in AOl where they find that extended synchrotron compo- 
nents very rarely contribute more than 1 mjy to the 850 ^im 
flux-density. However, AOl note that it is possible that two 
of the radio galaxy detections in our sample (3C 470 and 
6C 0820-1-3642) become non-detections after correcting for 
synchrotron emission. We do not attempt to make any cor- 
rections for this in this paper, due to the large uncertainties 
involved in extr apolating the spectrum from radio frequen- 
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cies. In Sec. 4.2 we will discuss this in terms of the angular 
extent of the radio sources with respect to the size of the 
SCUBA beam. Of course, radio quasars and radio galax- 
ies would suffer equally from lobe or hotspot emission, so 
it cannot explain the difference in submillimetre properties 
between the two types of object. 

The emission from radio cores and jets is beamed and 
Doppler boosting makes it highly orientation dependent. In 
the standard unification scheme, radio quasars have their 
jet axes closer to our line-of-sight and the enhanced Doppler 
boosting means that quasars have brighter cores than radio 
galaxies. Note that because the 3C and 6C radio samples 
were selected at low frequencies, where extended emission 
dominates, the quasars in our sample are lobe-dominated 
and not like the core-dominated blazars which are known 
to have fairly flat (0 ^5 a ^ 0.7) spectra extending into 
the submillimetre regime (e.g. Gear et al. 1994). There have 
been few investigations of the radio-millimetre spectral in- 
dices of lobe-dominated quasar cores. This is due largely 
to the fact that their 5 GHz flux-densities are typically in 
the range 10 — 100 mJy and, given that their core spectra 
are not flat, most have been below the sensitivity limit of 
many millimetre and submillimetre telescopes. In a study of 
six lobe-dominated quasars, Antonucci, Barvainis & AUoin 
(1990) flnd that the typical spectral index between 22 GHz 
and 230 GHz is a « 1.0. Van Bemmel & Bertoldi (2001) 
found the 5-250 GHz spectral indices of three radio quasars 
to all be « 1. Athreya et al. (1997) compare radio galaxy, 
lobe-dominated and core-dominated quasar cores and pre- 
dict that most lobe-dominated quasar cores are flat a « 
at GHz frequencies, but steepen at ~ 30 GHz (rest-frame). 
In Sec. 2 we showed that the core spectra of 3C 268.4 and 
3C 298 are also likely to have a ~ 1. Another one of our 
quasar sample, 3C 280.1, has an IRAM 1.25mm flux-density 
of 3.5 ± 1.7 (Andreani et al. 2001). This suggests its core 
spectral index from 5 GHz to 1.25 mm is q « 0.7, again in- 
dicating that the 850(i.m emission must contain a signiflcant 
contribution from dust. The one example of a flat {a ~ 0.3) 
millimetre spectral index for a lobe-dominated radio source 
is the unusual radio galaxy 6C 0902-1-34 (as discussed by 
AOl). The radio properties of this object suggest its jet axis 
is at a small angle to our line-of-sight and dust obscuration 
prevents us from observing the nuclear source and classi- 
fying it as a quasar (Carilli 1995). Therefore, the available 
evidence suggests that the core millimetre spectral indices of 
(lobe-dominated) radio quasars are mostly steep (« 1), but 
with a significant scatter such that some sources will have 
flatter core spectra. 
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Figure 4. 850 yxm flux-density Ssso against radio core flux- 
density at 5 GHz Score for radio quasars and radio galaxies (sym- 
bols as in Fig. 2). The curves correspond to the predicted con- 
tribution of core synchrotron radiation to Sgso as a function of 
Score and core spectral index a. 



In Fig. ^ we plot Sgso against the radio core flux- 
densities at 5 GHz given in Table 1. As expected, the radio 
quasars have cores which are brighter than those of the radio 
galaxies by a factor ~ 10. In Fig. ^ we also plot lines corre- 
sponding to the expected sub-mm flux-density of the core as 
a function of core flux-density for a range of spectral indices 
between 5 GHz and 850 yim. 3C 298 has a core flux a factor 
of ten greater than any of the other quasars in our study. 
As we saw in the previous section, its core almost certainly 
makes a significant contribution to its 850 \im fiux of « 10 
mJy. If the five quasars on the a = 0.5 line all have spectra 
with a ~ 1 then the synchrotron contribution to their 850 
[im fiuxes is ~ 1 mJy, well below their detected fluxes with 
the exception of 3C 205. Flatter spectral indices of a ~ 0.5 
would be necessary for the cores to be responsible for the 
detections of all these quasars. In reality, it is likely that 
there is a spread in the GHz to millimetre spectral indices 
of the quasar cores such that a few of these quasar submil- 
limetre detections might well be infiuenced by synchrotron 
emission. 

Five of the radio quasars (3C 191, 3C 268.4, 3C 298, 
3C 318 and 6C 1045+3513) are detected at shorter far- 
infrared wavelengths and from their SEDs we are confident 
that they are dusty hyperluminous objects (Sec. 2). For the 
other quasars, none of the 450 yun non-detections are sen- 
sitive enough for us to be able to exclude the presence of 
a thermal dust spectrum being responsible for the observed 
850 [im emission. 3C 181 and 3C 432 have very weak cores, 
which are steep (a > 0.5) at GHz frequencies and there- 
fore their cores cannot be signiflcant at 850 |i.m, so we con- 
sider these as further secure detections of dust. The IRAM 
1.25mm flux of 3C 280.1 suggests that its core component 
cannot account for all the observed 850 ^.m fiux. The remain- 
ing two submillimetre bright quasars (3C 205 and 3C 270.1) 
lie close to the a = 0.5 line in Fig. ^ and we cannot yet be 
certain that we have detected dust in these objects. Given 
the available evidence of the likely range of core spectra, 
we suspect that one or two of these submillimetre detec- 
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tions will be dominated by synchrotron emission. Sensitive 
millimetre observations of these sources are needed to dif- 
ferentiate between synchrotron and dust. We thus repeated 
the statistical tests performed in Section 2, now assuming, 
conservatively, that the submillimetre detections of both of 
these quasars are due solely to synchrotron emission: the 
detections were replaced with non-detections of 0.0 ± 1.0 
mjy. The mean quasar Ssso is now 5.2 ± 1.3, still a factor 
of 4 greater than the mean radio galaxy ^gso of 1.2 ± 0.6 
The survival analysis tests confirm that a statistical differ- 
ence between the submillimetre fluxes of radio galaxies and 
quasars remains at 95% significance. 



4 WHY DO SOME RADIO QUASARS HAVE 
HIGH SUBMILLIMETRE LUMINOSITIES? 

We have shown that core synchrotron emission may make a 
small contribution to the 850(i.m fluxes of 3C radio quasars, 
but cannot explain all the submillimetre difference between 
quasars and radio galaxies. In Fig. ^ we showed that there 
is no systematic difference in the radio luminosities of the 
two types of object. In this section we consider whether there 
are any other differences between radio galaxies and quasars 
which may help to explain the difference in their submillime- 
tre luminosities. The key to following our arguments is that 
we feel that in any realistic unifled scheme we must go be- 
yond the simplest model in which the probability of viewing 
a quasar nucleus is unity for a viewing angle 9 ^ Strans , with 
^'trails constant, and zero at 6 > 6'trans. 

4.1 The effect of the intrinsic luminosity of the 
quasar 

As discussed in Sec. 1 the idea that radio galaxies and ra- 
dio quasars with similar low-frequency radio luminosities 
should have similar distributions in a second orientation- 
independent quantity is not quite true. For example, Simp- 
son (1998) has quantifled the factor by which quasars 
should have higher [O III] line strengths than radio galaxies, 
matched in radio luminosity, in the context of the receding- 
torus model for AGN (Lawrence 1991). Considering a pop- 
ulation of quasars with a spread in optical luminosity, those 
quasars with higher optical luminosities have a larger torus 
opening angle due to the radius at which dust is sublimated 
being larger. Taking into account the scatter in the correla- 
tion between radio and optical luminosity, this means that 
sources with higher optical luminosities can be viewed within 
a larger opening angle and are thus more likely to be classi- 
fied as quasars. An additional effect is that the larger open- 
ing angles allow a greater fraction of the ionizing continuum 
to escape, however this is not the dominant effect in this 
case. Simpson (1998) showed that this leads to a sample of 
radio quasars having a mean optical luminosity, as probed 
by their [O III] line luminosity, approximately twice that of 
a sample of radio galaxies of the same radio luminosity. 

In Fig. 1^ we plot Ssso against the [O III] emission line 
luminosity -£/[oni], which is emitted isotropically, believed 
to be produced outside any obscuring torus, and which is 
a good indicator of the ionizing luminosity (Simpson 1998). 
Correlation tests show that Sgso and I/[oiii] for the radio 
quasars and radio galaxies are positively correlated at 95% 
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Figure 5. 850 flux-density 5g5o against narrow [OIII] emis- 
sion line luminosity i[oiii] for radio quasars and radio galaxies 
(symbols as in Fig. 2). 



significance. It is also apparent that the quasars do typi- 
cally have higher L[oiii] than the radio galaxies. For the 
quasars the mean log(I/[oiii]/W) is 36.9 ± 0.1 and for the 
radio galaxies it is 36.2 ± 0.2. Thus there is a factor of five 
difference in the mean [O III] luminosities of the two types 
of object in our sample. A similar difference between -£/[oiii] 
for the highest radio luminosity quasars and radio galaxies 
is also evident from Fig. 5 of Jackson & Rawlings (1997). 
In that study, all the objects have [O III] line observations 
made with the same wide slit, thereby eliminating uncer- 
tainties due to inhomogeneous datasets. This difference in 
emission line luminosity (and consequently ionizing lumi- 
nosity) could explain the difference of a factor a few in the 
average luminosities of quasars and radio galaxies in the 
mid- and far-infrared in studies with IRAS and ISO (Heck- 
man et al. 1994; Hes et al. 1995; van Bemmel, Barthel & de 
Graauw 2000; Meisenheimer et al. 2001), since this radia- 
tion is certainly dominated by nuclear emission reprocessed 
by the dust torus. 

Could such an effect also be responsible for the similar 
factor difference in the submillimetre luminosities of radio 
quasars and radio galaxies? This would require a tight cor- 
relation between the submillimetre and optical continuum 
emission from quasars. The analysis of Simpson (1998), pre- 
dicting a factor of two difference in the [O III] luminosities 
of radio quasars and radio galaxies, used the known scatter 
of 0.6 dex in the radio-optical correlation (Serjeant et al. 
1998; Willott 2000). Therefore, we would expect the scat- 
ter exhibited in the radio-submillimetre correlation to be of 
the same order if the submillimetre and optical continua are 
tightly correlated. A power law fit to the marginal correla- 
tion between L151 and Sgso in Fig. 2 gives a linear slope with 
a scatter of 0.3 dex. The scatter may be slightly larger than 
this because of the large number (40%) of non-detections (in 
the above analysis non-detections were set to their la upper 
limits). However, within the limited range in L151 available 
in our study, we find the scatter in the radio-submillimetre 
correlation to be no larger than that of the radio-optical 
correlation. 
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A tight correlation between submillimetre and optical 
luminosities could probably only come about by direct heat- 
ing of the dust responsible for the submillimetre emission by 
the AGN (as is known to be the case for the mid-infrared 
emission, and has been discussed previously by Willott et 
al. 2000b). In this case the scatter could be as small as that 
seen in the correlation between optical and mid-infrared lu- 
minosities of optically-bright, low-redshift quasars, i.e. 0.2 
dex from the study of Rowan- Robinson (1995). If the far- 
infrared emission is due to dust heated by an intense star- 
burst rather than by the quasar, a larger scatter is predicted; 
indeed the larger scatter, of 0.4 dex, in the far-infrared ver- 
sus optical luminosity plot from Rowan- Robinson (1995) was 
used as an argument that: (i) the far-infrared emission from 
the low-redshift quasars is dominated by starburst-heated 
dust; and (ii) the scaling relation between optical and far- 
infrared luminosities arises from both quantities being corre- 
lated with a third variable, i.e. the black hole and/or galaxy 
mass. However, in the quasar-heating hypothesis, the far- 
infrared luminosity is strongly dependent upon the extended 
dust distribution (e.g. Granato & Danese 1994) and would 
indeed be expected to show a larger scatter than the mid- 
infrared-optical ratio. It may seem surprising that objects 
with such large quantities of dust have such a small scatter 
in these ratios because of the effects of dust-reddening of 
the optical emission. However, in both the Rowan- Robinson 
sample and the radio-loud sample discussed in this paper, 
objects would not be classified as quasars if they were heav- 
ily reddened {Ay ^ 3). 

Finally in this section, we note that despite the quan- 
titative agreement between the predictions of the receding- 
torus model and the submillimetre difference between radio 
quasars and radio galaxies, other physical effects seem likely 
to be important. Note, for example, that the high submil- 
limetre luminosity of 6C 1045-f3513 is an anomaly in the 
receding-torus model, since it has a value of I/[oiii] below 
that of all the 3C radio galaxies. 

4.2 The effect of the radio source age 

Fig. ^ plots 5*850 against the projected linear sizes D of 
the radio sources. According to unified schemes, quasars 
will typically have projected linear sizes a factor of « 1.6 
smaller than those of radio galaxies, due to the smaller an- 
gles between their jet axes and our line-of-sight. (assuming 
^trans = 53° from Willott et al. 2000a). Taking account of 
this effect, the distributions in D for the quasars and ra- 
dio galaxies are fairly similar. However, a striking feature 
of this plot is that the quasars appear to show an anti- 
correlation between Sgso and D. The correlation tests per- 
formed (Kendal's Tau and Cox's) show these quantities to 
be anti-correlated with significances of 95 and 99%, respec- 
tively. Including the radio galaxies as well, and accounting 
for projection effects, the significance of this correlation is 
only marginal (92% and 50%), largely because of the intro- 
duction of a large number of submillimetre upper limits. 

An important point to consider here is whether this 
anti-correlation could be due to contamination of the 
sub-millimetre emission by synchrotron from the lobes or 
hotspots. This is because the SCUBA beam at 850 [i-m has a 
13 arcsec FWHM, corresponding to a physical size at these 
redshifts of 110 kpc. Hence sources with D less than this 
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Figure 6. 850 flux-density Sgso against radio source pro- 
jected linear size D for radio quasars and radio galaxies (symbols 
as in Fig. 2). There is a significant anti-correlation present for 
the quasars. The upper axis shows how the ages of radio sources 
approximately relate to D, according to the model of Willott et 
al. (1999) (see Sec. 4.2 for details). 

value will have their hotspots within the SCUBA beam. 
However, one should note that the five radio quasars with 
detections at shorter wavelengths indicated by the crosses 
are also mostly the smallest quasars. In these objects we 
have evidence for a thermal dust spectrum and high far- 
IR luminosities. It can also be seen that five radio galaxies 
have D < 100 kpc and these have significantly smaller sub- 
millimetre fiuxes than the quasars. In the orientation-based 
unified scheme there is no reason why quasars should have 
brighter extended synchrotron emission than radio galaxies. 
Therefore, although we cannot rule out some lobe/hotspot 
contamination in the smallest sources it cannot be the driv- 
ing force behind the observed correlation. 

3C 298, 3C 318, 6C 1045-^3513, and more marginally 
3C 191, are all compact steep-spectrum (CSS) radio sources, 
which are believed to be small due to their jets being trig- 
gered relatively recently (~ 10^ years; Fanti et al. 1995). 
Previous studies of the far-infrared properties of CSS sources 
have shown them to be indistinguishable from larger, pre- 
sumably older, radio sources (Heckman et al. 1994; Fanti et 
al. 2000). However our study is quite different from these 
which were limited to low redshifts (and hence luminosi- 
ties) due to the sensitivity limitations of IRAS and ISO, and 
it is plausible that the differences only occur at higher lu- 
minosities, higher redshifts and/or longer wavelengths (e.g. 
because the emission is dominated by starburst-heated cold- 
dust emission rather than AGN-heated hot-dust torus emis- 
sion) . 

It is also interesting to note that 3C 298, 3C 318 and 
6C 1045+3513 are reddened quasars, indicating significant 
quantities of dust in their host galaxies. Also, 3C 191, 
3C 268.4 and 3C 298, show very strong CIV associated 
absorption and 6C 1045-1-3513 has a CIV broad absorp- 
tion system. Baker et al. (2001) show that there is an 
anti-correlation between the strength of associated absorp- 
tion systems and radio linear size and also that the redder 
quasars tend to have the strongest line absorption. It has 



A submillimetre difference between radio galaxies and radio quasars 9 



also been noticed that, amongst radio-loud quasars, broad 
absorption lines are more commonly observed in the smallest 
radio sources (Becker et al. 2000). These observations indi- 
cate an evolutionary picture where the absorbing material 
is cleared as the radio sources get older (Baker et al. 2001). 
This is similar to the evolutionary scheme of Sanders et al. 
(1988) where quasars begin their lives in a dust-enshrouded 
phase. 

The projected linear size D of a radio quasar or radio 
galaxy can be related approximately to the time since the 
jets were triggered t, since the relationship between these 
variables is only weakly dependent upon the radio lumi- 
nosity, and the objects under consideration span only 1.5 
dex in this quantity. Using the simple model of radio source 
evolution in Willott et al. (1999) we can, at least approxi- 
mately, transform the linear size axis in Fig. ^ to the time 
since the quasars were triggered (the radio galaxies have 
ages on average a factor of 1.6 smaller at a given D than 
shown here, due to projection effects). Fig. ^, and the anti- 
correlation it exhibits, shows that the most submillimetre- 
luminous quasars had their radio sources triggered < 10^ 
years ago. This means that the processes controlling the 
submillimetre/far-infrared luminosity Lfir are synchronized 
with the jet-triggering event and that Lpir declines on a 
~ 10^ yr timescale. 

Although high values of Lpir for an AGN are generally 
interpreted as evidence for high rates of star formation, it 
is possible that, as explored in Sec. 4.1 (see also Willott et 
al. 2000b; Rowan- Robinson 2000), even the cool dust is di- 
rectly heated by the AGN. Any evidence for star formation 
is indirect, being based on the large gas masses observed in 
CO in a handful of cases, but mostly deduced from the large 
dust masses required by any submillimetre-detected thermal 
source (e.g. Hughes, Dunlop, & Rawlings 1997). With a stan- 
dard dusty disc model (e.g. Efstathiou & Rowan-Robinson 
1995) there are difficulties in heating the cool dust at large 
distances from the nucleus, due to the fact the disc is op- 
tically thick and shields the distant dust from the ionizing 
radiation. A more efhcient way of heating large amounts of 
cool dust by an AGN, come from models where the dust has 
a more isotropic distribution with a very high covering factor 
(e.g. Granato, Danese & Franceschini 1996). This may oc- 
cur preferentially in recently-triggered AGN if the dusty gas 
fueling the nucleus is not in a regular torus distribution. In 
addition, models for shocks caused by advancing radio com- 
ponents suggest that dust grains cannot survive for much 
longer than ^ lO^yr in such environments (De Young 1998). 
Jarvis et al. (2001b) summarise other evidence that expand- 
ing radio sources modify the physical conditions of their 
gaseous environments as the radio-emitting components ex- 
pand outward after jet triggering. There are therefore two 
reasons why in a scenario in which quasar-heated dust dom- 
inates, one expects Lfir to decline on a ~ lO^yr timescale. 
Note that the fact that the brightest submillimetre sources 
are reddened quasars with line absorption is consistent with 
this picture, since with this dust distribution there is a high 
probability of our line-of-sight to the nucleus intercepting a 
dust cloud. It is of course possible that the dust distribution 
could be influenced by previous episodes of star- formation. 

It is also possible that the submillimetre emission in 
the hyperluminous sources is dominated by starburst-heated 
dust, in which case the implication is that the jet-triggering 



process is accurately synchronized with an intense burst of 
star formation. The timescale for such star formation events 
would have to be ^ 10^ years for the far-infrared (300 \xm) 
luminosity to peak « 10^ yr after the onset of the star- 
burst (Efstathiou, Rowan- Robinson & Siebenmorgen 2000). 
Line ratios from mid-infrared spectroscopy of nearby star- 
burst galaxies suggest that the duration of the bursts in such 
galaxies are ^ 10'^ yr (Thornley et al. 2000) and for ultralu- 
minous infrared galaxies ~ 2 X 10^ yr (Genzel et al. 1998). 
Thus the starburst timescales are quite consistent with the 
hypothesis that a starburst and a jet-producing AGN were 
triggered synchronously, and the anti-correlation seen in Fig. 

i 

A starburst of such a brief duration, forming stars at 
^ 1000 Mq yr~^ as implied by the submillimetre luminos- 
ity, would form ~ 10^" Mq of stars over its lifetime. The 
total stellar mass of a powerful radio galaxy with a lumi- 
nosity of 5L* (Jarvis et al. 2001a) is 5 x 10"Mq. Therefore 
such a star formation episode would form only a few per 
cent of the total stellar mass in a massive elliptical galaxy 
and therefore this phase is not a major event in the star- 
formation histories of the quasar host galaxy. This is con- 
sistent with other evidence, such as the small scatter in the 
K — z relation (Jarvis et al. 2001a) and the existence of very 
red radio galaxies at 1 < z < 2 (Dunlop et al. 1996; Dunlop 
1999; Willott, Rawlings & Blundell 2001), that the bulk of 
the stars in powerful radio galaxies were formed at higher 
redshifts of 2 > 3. 

Regardless of whether the dust is heated by a quasar or 
a starburst or a mixture of the two, it seems that Fig. ^ pro- 
vides some support for the Blundell & Rawlings (1999) con- 
jecture that radio source youth may be an important factor 
in explaining the high submillimetre luminosities of distant 
radio galaxies. An important question is how this finding 
relates to the properties of the more distant (2 > 2) systems 
studied by AOl. We again defer a more detailed discussion to 
a future paper (Rawlings et al. in prep.) but note here that 
the fact that the submillimetre luminosities of the 2 < 2 ra- 
dio quasars are similar to the radio galaxies at 2 > 2 is not, 
by itself, a strong argument that their host galaxies are in 
similar evolutionary states. It is certainly plausible that the 
higher-redshift objects are forming significant fractions of 
the stellar mass of a giant elliptical galaxy, requiring simply 
that the duration of the starburst is ';S> lO'^ yr. Direct mea- 
surements of gas mass, e.g. by CO observations, would be 
very important for determining the evolutionary states and 
the main heating source of the cool dust for the 2 ~ 1.5 ob- 
jects studied here. It is hard with submillimetre continuum 
measurements alone to determine the relationship between 
these objects and the high-redshift radio galaxies studied 
by AOl. The discovery of huge reservoirs of molecular gas 
in some very high-redshift radio galaxies (Papadopoulos et 
al. 2000; Alloin, Barvainis & Guilloteau 2000) supports a 
starburst-heating model for the bulk of the submillimetre 
continuum in these objects, and the lack of such reservoirs 
at 2 ~ 1.5, if observed, would suggest that quasar-heating 
might become more prevalent at lower redshifts. 
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5 CONCLUDING REMARKS 
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We have performed a submillimetre survey of radio quasars 
in the redshift range 1.37 < z < 2 and detect ten out 
of eleven quasars, implying submillimetre/far- infrared lu- 
minosities ^4 X IO^^Lq. We argue that only one or two 
of these detections are dominated by core synchrotron emis- 
sion. For five quasars, detections at shorter wavelengths con- 
firm the presence of a thermal continuum and a large mass 
(~ 10*Mq) of cool dust. 

We find the radio quasars have higher submillimetre 
luminosities by a factor of ~ 5 than radio galaxies of the 
same radio luminosity and redshift, a factor which cannot be 
reduced below ~ 2 by accounting for possible synchrotron 
contamination. At 2: ~ 1.5 radio quasars are intrinsically 
brighter than radio galaxies at submillimetre wavelengths. 
This argues against the simplest unified scheme for radio 
sources in which the probability of viewing a quasar nucleus 
is unity for a viewing angle 6 ^ 6'trans, with 6trans constant, 
and zero &t 9 > Strans. It is, however, in quantitative agree- 
ment with a receding-torus model in which ^trans depends 
on quasar optical luminosity, but only if there is a close rela- 
tionship between optical luminosity and Lfir. The implica- 
tion is that quasar-heated dust dominates Lfir for powerful 
radio quasars at z ~ 1.5. 

There is a significant anti-correlation between Lfir and 
radio source age in the sense that hyperluminous quasars 
tend to be associated with young (< 10^ yr) sources at 
z ~ 1.5. This means that the processes controlling Lfir are 
synchronized with the jet-triggering event and that Lfir de- 
clines on a ~ 10^ yr timescale. There is evidence for a small 
amount of obscuration in the most submillimetre-luminous 
quasars from reddening (3C 298, 3C 318 and 6C 10454-3513) 
and/or prominent associated absorption (3C 191, 3C 268.4, 
3C 298) or broad absorption lines (6C 1045-1-3513). We spec- 
ulate that shortly after jet triggering, dust is more widely 
distributed than at later times. This distributed dust could 
be heated by the quasar nucleus and/or a short-lived syn- 
chronized starburst and may be destroyed by shocks asso- 
ciated with the expanding radio source. Any such starburst 
would form only a few per cent of the total stellar mass in 
agreement with the evidence for dominant old stellar popu- 
lations in 2 ~ 1.5 radio galaxies. 
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APPENDIX A: 6C 1045+3513 - A 
SUBMILLIMETRE-LUMINOUS REDDENED 
BAL QUASAR 

In this appendix we present a new optical spectrum of 
6C 1045+3513, which is one of the brightest radio quasars 
at submilhmetre wavelengths found in this paper (5*850 = 
9.18+1.10 mjy). It is also detected at the Sa level with the 
short wavelength array (5*450 = 21.3 ± 6.8 mJy). 

6C 1045+3513 was observed with the ISIS double-beam 
spectrograph on the William Herschel Telescope on April 16 
2001. Conditions were photometric with 1.0 arcsec seeing 
and a 1.0 arcsec slit was used giving a dispersion resolution 
of 8 A. The telescope was pointed at the position of the com- 
pact radio source (10:48:34.25 +34:57:25.0 - J2000.0). The 
total integration time of the observation was 1200 seconds. 
The data were reduced in the IRAF package using the same 
method as in Willott et al. (1998). The red and blue arm 
spectra were joined together, averaging over 35 A where 
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Wavelength (A) 

Figure Al. Optical spectrum of the reddened quasar 
6C 1045+3513 with emission lines labelled. The CIV broad ab- 
sorption feature is also marked. The absorption features marked 
with © are due to atmospheric absorption and that marked 'di' 
is an artefact caused by the dichroic. 



the signal-to-noise ratio of the spectra were approximately 
equal (6200 A). 

The target is clearly identified in the two dimensional 
spe ctru m and appears to be spatially unresolved. In Fig- 



Al 



we show the optical spectrum extracted in a 1 arc- 
sec aperture. The object has a red continuum and several 
strong emission lines (see Table Al). The red continuum 
has a spectral index of 2.3 and is well fit by a quasar spec- 
trum reddened by Ay = 1.5 (assuming a Milky-Way ex- 
tinction curve). The positions of the emission lines are con- 
sistent with a redshift for 6C 1045+3513 oi z = 1.604, ex- 
tremely close to the tentative value reported by Rawlings et 
al. (2001). The CIV and Hell emission lines are both narrow 
(FWHM f» 2000 km s"^). The Mg II and C HI] emission lines 
appear to have narrow central profiles with broad wings, 
which together with the weak or absent broad wings in C IV 
and He II, are consistent with broad-line region reddening of 
the same amount as the continuum fitting suggests. 

There is very high equivalent width C IV absorption on 
the blue wing of the emission line. This feature is very broad 
with a width of at least 6000 km s~^. There is no corre- 
sponding broad absorption of Mgll, so 6C 1045+3513 is 
classified as a high-ionization broad absorption line (BAL) 
quasar. Unfortunately, our data are not sensitive enough or 
with sufficient spectral resolution to determine the structure 
of the absorption trough. 6C 1045+3513 is one of the most 
radio-luminous quasars to show the BAL phenomenon, with 
a radio luminosity similar to that of the BAL QSO FIRST 
J101614.3+5209I6 (Gregg et al. 2000) . 

6C 1045+3513 is thus remarkably similar to 3C 318 
(Willott et al. 2000b). Both objects have a small degree 
of reddening which suppresses the rest-frame UV broad 
line strengths to below those of the narrow lines. Both are 
also very luminous submilhmetre sources (with detections at 
both 850 (i.m and 450 \im) and both are compact {D < 10 
kpc) radio sources with similar estimated source ages. A 
possible reason for these similarities is discussed in Sec. 4.2. 
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Line 


(A) 




FWHM 

(km s^-^) 


Flux/10~i^ 
(W m-2) 


CIV A1549 


4029 ± 4 


1.601 


2000 


18 ± 4 


He II A1640 


4263 ± 6 


1.599 


2400 


3.6 ± 1.3 


CIII] A1909 


4906 ± 6 


1.597 


4300 


8.2 ± 1.0 


Mgll A2799 


7291 ± 2 


1.605 


2900 


9 ± 4 


[NeV] A3426 


8920 ± 3 


1.604 


1300 


4.6 ± 1.6 



Table Al. Emission line data for 6C 1045+3513 from our optical 
spectrum. Fluxes have been aperture corrected. 



